Observation of anisotropic diamagnetism above the superconducting transition in 
iron-pnictide Ba 1 _ ;c IC r Fe2As 2 single crystals due to thermodynamic fluctuations 

J. MosqueiraQ JD. Dancausa, and F. Vidal 
LBTS, Facultade de Fisica, Universidade de Santiago de Compostela, E-15782 Santiago de Compostela, Spain 

S. Salem-Sugui, Jr. 

Institute de Fisica, Universidade Federal do Rio de Janeiro, 21941-972 Rio de Janeiro, RJ, Brazil 



o 

(N 



in 

(N 



P3 
O 
O 

u 

Q-i! 



o 

> 

OO 

in 

d 
o 



X 



A.D. Alvarenga 

Instituto Nacional de Metrologia Normalizacao e Qualidade Industrial, 25250-020 Duque de Caxias, RJ, Brazil 

H.-Q. Luo, Z.-S. Wang, and H.-H. Wen 

National Laboratory for Superconductivity, Institute of Physics and national Laboratory for Condensed Matter Physics, 

P.O. Box 603, Beijing 100190, People's Republic of China 
(Dated: March 28, 2011) 

High resolution magnetization measurements performed in a high quality Bai_ a: K I Fe2As2 single 
crystal allowed to determine the diamagnetism induced above the superconducting transition by 
thermally activated Cooper pairs. These data, obtained with magnetic fields applied along and 
transverse to the crystal ab layers, demonstrate experimentally that the superconducting transition 
of iron pnictides may be explained at a phenomenological level in terms of the Gaussian Ginzburg- 
Landau approach for three-dimensional anisotropic superconductors. 
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I. INTRODUCTION 



The recent discovery of superconductivity in iron- 
pnictides^ has generated an intense research activity, in 
part because they provide an unexpected and very in- 
teresting scenario to study the interplay between mag- 
netism and superconductivity— The interest for these 
compounds is enhanced by the fact that they share many 
properties with the high-T c cuprates (HTSC), as the lay- 
ered structure, a similar evolution of the superconducting 
parameters with doping, and the proximity to a magnetic 
transition^ As it is still the case of the HTSC, the pair- 
ing mechanism in superconducting iron pnictides is not 
yet known, making the phenomenological descriptions of 
their superconducting transition a central issue of the 
present researches on these compounds. As first estab- 
lished in the pioneering experiments by M. Tinkham and 
coworkers in low-T c metallic superconductors^ a power- 
ful tool to probe at a phenomenological level the own 
nature of a superconducting transition is the diamag- 
netism induced above T c by thermally activated Cooper 
pairs^ The fluctuation diamagnetism above T c (some- 
times called precursor diamagnetism) was also early used 
to characterize at a phenomenological level the super- 
conducting transition in both optimally-doped^ and 
underdoped 7 high-T c cuprates. However, in spite of its 
interest the precursor diamagnetism received little atten- 
tion due to both the relative smallness of the supercon- 
ducting fluctuation effects above T c in these compounds 
and, mainly, to the need of high quality single crystals, 
with a relatively sharp superconducting transition. 8 



ducting iron-pnictides of the fluctuation-induced dia- 
magnetism above T c , which were possible by applying 
high resolution magnetometry to a Bai_ a: K a ;Fe2As2 sin- 
gle crystal with an excellent chemical and structural 
quality. These data allow an experimental demonstra- 
tion that the superconducting transition of these ma- 
terials may be explained at a phenomenological level 
in terms of the Gaussian Ginzburg-Landau approach 
for three-dimensional anisotropic superconductors (3D- 
AGL). This conclusion excludes phase incoherent super- 
conductivity above T c in iron-pnictides, a long standing 
but still at present debated issue of the HTSC physics i^iIS 
Also, it provides a direct check of the fluctuations dimen- 
sionality: while bulk low-T c superconductors (LTSC) be- 
have as three-dimensional (3D) and most HTSC as two- 
dimensional (2D)p^i it has been recently proposed that 
iron-pnictides are in an intermediate regime as a conse- 
quence of their moderate anisotropyjii In relation with 
this, another interesting aspect of our work is that the 
fluctuation diamagnetism is probed with magnetic fields 
applied along the two main crystallographic directions. 
This allows to accede directly to the superconducting 
anisotropy factor, but also to study an interesting issue 
recently observed in other moderately anisotropic super- 
conductor (NbSe2), as it is the anisotropy of non-local 
electrodynamic effects^ 



In this paper we report measurements in supercon- 
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II. EXPERIMENTAL DETAILS AND RESULTS 
A. Sample fabrication and characterization 

The Bai_ :c K 2 ,Fe2As2 sample (with x — 0.28) used in 
this work is a ~ 1.1x0.7x0.1 mm 3 (430 /zg) single crystal. 
Details of its growth procedure and characterization may 
be seen in Ref. [H. We have chosen such a small crystal 
because it presents a well defined transition temperature 
which is essential to study critical phenomena. The mea- 
surements were performed with a magnetometer based 
in the superconducting quantum interference (Quantum 
Design, model MPMS-XL). To position the crystal with 
the magnetic field, H, applied parallel to the crystal ab 
layers, we used a quartz sample holder (0.3 cm in diame- 
ter, 22 cm in length) to which the crystal was glued with 
GE varnish. Two plastic rods at the sample holder ends 
(~ 0.3 mm smaller than the sample space diameter) en- 
sured that its alignment was better than 0.1°. For the 
measurements with H _L ab we made a groove (~ 0.3 
mm wide) in the sample holder into which the crystal 
was glued also with GE varnish. The crystal alignment 
was checked by optical microscopy to be better than 5°. 
This allowed to determine the anisotropy factor from the 
anisotropy of the precursor diamagnetism with a ~ 0.5% 
uncertainty^ 

As a first magnetic characterization, in Fig. 1 it is 
presented the temperature dependence of the zero-field- 
cooled ZFC and field-cooled FC magnetic susceptibilities, 
measured with a magnetic field of 0.5 mT perpendicular 
to the ab planes. These measurements were corrected 
for demagnetizing effects by using as demagnetizing fac- 
tor D = 0.78, which was obtained by approximating the 




FIG. 1. (Color online) Temperature dependence of the FC 
and ZFC magnetic susceptibility (already corrected for de- 
magnetizing effects) obtained with a 0.5 mT magnetic field 
perpendicular to the ab layers. T c and its uncertainty were 
estimated from the maximum of d X ZFC (T)/dT (solid line), 
and from the high-temperature half-width at half-maximum 
(this eludes the possible widening associated to the transit 
through the mixed state). 



crystal shape by the inscribed ellipsoid. A proof of the 
adequacy of this procedure is that the ZFC magnetic sus- 
ceptibility in the Meissner region is close to the ideal 
value (-1) well below T c . From these curves we estimated 
T c ss 33.2 ± 0.2 K, which attests the excellent quality of 
the sample, and allows to study the fluctuation-induced 
magnetization in almost the whole temperature range 
above T c . Such a small transition width is in good agree- 
ment with the one determined from the resistive tran- 
sition of crystals with the same composition and grown 
following the same procedure.— 



B. Measurements of the fluctuation diamagnetism 
around T c and background subtraction 

To measure the effect of superconducting fluctuations 
above T c (which according to the 3D-AGL approach, see 
below, is about 10~ 6 -10 -7 emu), we used the reciprocat- 
ing sample option. It produces sinusoidal oscillations of 
the sample about the center of the detection system and 
improves the resolution by about two orders of magnitude 
with respect to the conventional DC option. Data were 
taken by waiting ~ 3 min for complete stabilization after 
the temperature was within 0.5% the target temperature. 
For each temperature we averaged six measurements con- 
sisting of 15 cycles at 1 Hz frequency. The final resolution 
in magnetic moment, m, was in the ~ 10~ 8 emu range. 

An example of the as measured m(T) data above T c 
is presented in Fig. 2(a). As may be clearly seen al- 
ready in this figure, the m(T) curves present a rounding 
just above T c (a factor ~ 3 larger in amplitude when 
H _L ab), which extends few degrees above T c . In view of 
the sharp low-field diamagnetic transition (Fig. 1) such 
a rounding cannot be attributed to T c inhomogeneities, 
and is a first evidence of the presence of an anisotropic 
precursor diamagnetism in these materials. In the detail 
around T c presented in Figs. 2(b) and (c), it is clearly 
seen that the reversible region extends few degrees be- 
low T c , allowing to study the critical fluctuation regime, 
where m(T) presents a similar anisotropy. 

The fluctuation-induced contribution was determined 
from the as measured m(T) by subtracting the back- 
ground magnetic moment, mg(T), coming from the crys- 
tal normal-state and in a much lesser extent from the 
sample holder (this last contribution was estimated to 
be in the 10~ 6 emu range). Such a background contribu- 
tion was determined by fitting the function 



ms(T) = a + bT - 



T 



(1) 



to the raw data in the temperature interval between 37 
and 50 K (a, b and c are free parameters). The lower 
bound corresponds to a temperature above which the 
fluctuation-induced magnetic moment is expected to be 
smaller than the experimental uncertainty (see below). It 
is worth noting that the use of other plausible functions 
for mg(T), for instance ms(T) = a + b/(T — c), leads to 
imperceptible changes in the scale of Fig. 2(a). 
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FIG. 2. (Color online) (a) Example of the temperature depen- 
dence up to 50 K of the as-measured magnetic moment for a 
magnetic field applied in the two main crystallographic direc- 
tions. The normal-state backgrounds (lines) were determined 
by fitting Eq. (QJ in the region indicated. The shaded ar- 
eas represent the precursor diamagnetism. (b) and (c) Detail 
around the superconducting transition. Solid (open) symbols 
were obtained under ZFC (FC) conditions. The diamagnetism 
above T c is unobservable in this scale. 



III. DATA ANALYSIS 

A. Theoretical background 

In view of the moderate anisotropy observed, the 
results of Fig. [5] were analyzed in the framework of 
the 3D-AGL approach. This theory predicts that the 
fluctuation-induced magnetization for magnetic fields ap- 
plied perpendicular or parallel to the crystal ab layers 
(Mj_ or M|| , respectively) may be obtained from the re- 
sult for 3D isotropic materials, Mj so , through^— 



M ± (H) = jM iso (H), 



and 
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where 7 is the superconducting anisotropy factor. This 
transformation was introduced by Klemm and Clemi^ 
and generalized by Blatter— and by Hao and Clemi£ to 
different observables and different regions in the H — T 
phase diagram. In particular, it was shown to be also 
valid for the fluctuation region above T C J£ For isotropic 
materials, the fluctuation magnetization above T c was 
calculated in Ref. in the framework of a Gaussian GL 
approach including a total-energy cutoff in the fluctua- 
tion spectrum^ By combining that result with Eqs. © 
and © it follows 



Mi = - 



k B Tfj, o H^ ab (0) 



arctan y/(c— e)/e 



arctan \J[c~ e)/c 



and 



My = 



M ± 
— 



(4) 



(5) 



where e — ln(T/T c ) is the reduced temperature, £ a b(0) 
the in-plane coherence length amplitude, c 0.55 the 
cutoff constant, 19 k B the Boltzmann constant, /j,q the 
vacuum magnetic permeability, and <po the magnetic flux 
quantum. These expressions are valid in the low mag- 
netic field limit, i.e., for H <C eH^O) and, respec- 
tively, H < e-ffjko) = £7^(0), where H^(0) = 
<^o/27r/io£ab(0) is the upper critical field for H _L ab ex- 
trapolated linearly to T = K. In absence of any cutoff 
(i.e., when c — > 00) Eq. (4) simplifies to 



M, = - 



Trk B T[i o H-/£ ab (0) _ 1 



/2 



(6) 



which corresponds to the well known Schmidt result^ 

Equations dU) and J5|) predict the vanishing of the pre- 
cursor diamagnetism at e — c (i.e., T ss 1.7T C « 56 K 
in the present case). This has been experimentally con- 
firmed in a number of HTSC and LTSCji2 but in the 
present case the small crystal size does not allow a quan- 
titative check of this onset temperature^ The difference 
between Eq. ^ and the conventional Schmidt approach 
Eq. ^ is dramatic at high reduced temperatures (typi- 
cally above e = 0.1). However, it is still as high as 15% 
for e as low as 0.01, which justifies the use of the total 
energy approach also at low reduced temperatures. 



B. Fluctuation diamagnetism in the Gaussian 
region above T c 

A detail of the fluctuation magnetization above T c (al- 
ready corrected for the normal state contribution and 
normalized by the applied field) is presented in Fig. 3 
for both H _L ab and H \\ ab. In this representation 
we have only included magnetic fields up to 2 T, which 
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are expected to be in the low magnetic field regime ex- 
cept in a narrow temperature interval close to T c . As 
expected, both M±/H and Mn/H are independent of 
the applied magnetic field except very close to T c , where 
finite-field effects are already observable The best 
fit of Eqs. ((4| and (O to the data in the low-field region 
is in good agreement with the temperature dependences 
and amplitudes observed. The values obtained for the 
two free parameters, £ a fc(0) = 1.5 nm (which leads to 
HodH^/dT = -4.2 T/K near T c ) and 7 = 1.8, are close 
to the ones in the literature^ These values lead to a 
transverse coherence length £ c (0) = £a&(0)/7 = 0.83 nm, 
larger than the Fe-layers periodicity length (^0.66 nm), 13 
which justifies the applicability of the 3D-AGL approach. 
The £ab(0) value is well within the ones in HTSC, and 7 
is just a few times smaller than in YBa2Cu307_,5. These 
similarities contrast with the differences between their 
precursor diamagnetism (allegedly unconventional in the 
HTSC)£ Finally, £ a t>(0) is much shorter than in the mod- 
erately anisotropic 2H-NbSc2, which could explain the 
non observation of non-local electrodynamic effects when 
H _L ab in the present caseJ^ 



Fluctuation diamagnetism in the critical region 
around T c 



As a check of consistency of the above analysis we 
studied the data in the critical region, closer to the in- 
dependent critical temperature, where the fluctuations 
amplitude is so large that the Gaussian approximation 
breaks down.— This region is bounded by the so-called 
Ginzburg criterion^ which according to the transforma- 
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FIG. 3. (Color online) Detail around T c of the temperature 
dependence of the fluctuation magnetization (over H) for sev- 
eral magnetic fields applied in the two main crystallographic 
directions. The lines correspond to the Gaussian 3D-AGL ap- 
proach in the low magnetic field limit (solid blue for H _L ab 
and dashed red for H \\ ab). For a better comparison, both 
lines are presented in (a) and (b). 
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FIG. 4. (Color online) 3D-GL scaling of the magnetization in 
the critical region for H _L ab (a) and H \\ ab (b). The solid 
blue line is the experimental scaling function for H _L ab, while 
the dashed red one corresponds to H \\ ab and was calculated 
from /j_ through Eq. ([lip, (c) and (d) H-T phase diagrams for 
H _L ab and, respectively, H \\ ab. The corresponding H C 2(T) 
lines were obtained from the superconducting parameters ob- 
tained in the analysis. The symbols are the lower bound of 
the reversible region, and the dashed lines the limits of the 
critical-region according to Eqs. Q and ([8]). 



tion to anisotropic materials may be written as 



T kT^{H)±T c 
for H _L ab, and 

T « TJ 1 (H) ± T r 



Ac£ c (O)0 o 
AnkBtJ-aH 



2/3 



2/3 



(7) 



(8) 



c [Ac£ c (0)7^>o 

for H || ab, where Ac is the specific-heat jump at T c . In 
this region, the 3D-GL approach in the lowest-Landau- 
level approximation predicts that M± (T, H) follows a 
scaling behavior, m± = f±(t±), the scaling variables 
being 2 ^ 



and 



" (iiT) 2 / 3 
T-T^{H) 



(10) 



^~ (HT) 2 / 3 ' 

where T^{H) = T c [l - H/H^(0)}. According to Eq. ©, 
Mil (T, H) should scale as 



m ll =/||( f ll) 



/i(*||7 2/3 ) 



v 5/3 
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the scaling variables being now 



Mi, 



(HT)V 3 



(12) 



to be adequate when H _L abr- but here its validity is 
extended to the case in which H II ab. 



IV. CONCLUSIONS 



and 



T - rj (g) 
(JTT) 2 /3 ' 



(13) 



where t}{H) = T c [l-H/H l } 2 (0)}, and ffj 2 (0) = 7#ci(0). 

In Figs. SJa) and (b) it is presented the scaling of the 
Mj_(T, H) and My (T, ff) data according to Eqs. (|9l>. (fTO)) . 
(H2) and ([13"]) . For that, t± and tii were evaluated by us- 
ing the H^(0) and 7 values obtained in the analysis of 
the precursor diamagnetism in the Gaussian region, and 
T c = 33.1 K (in good agreement with the T c value result- 
ing from Fig. 1). As may be clearly seen, both scalings 
are excellent and the relation between the observed scal- 
ing functions is in good agreement with Eq. (1111) . The 
applicability of the above scalings extends down to the 
lower bound of the reversible region. This is in agreement 
with Eqs. ([7]) and <j8j) , which when evaluated by using 
Ac/T c = 0.1 J/molK 2 (see Ref. [H and the supercon- 
ducting parameters resulting from the previous analysis, 
are in an unexpected agreement with the observed lower 
bound of the reversible region for both H directions [see 
Figs. 4(c) and (d)]. The 3D scaling was previously found 



The present experimental results and analysis repre- 
sent a solid evidence that the precursor diamagnetism 
and the magnetization in the critical region around T c 
in Bai_ a: K 2 ;Fe2As2 may be consistently explained at a 
phenomenological level in terms of a conventional GL 
approach for 3D anisotropic superconductors. This pro- 
vides a strong constraint for any microscopic theory for 
the superconductivity in iron pnictides, and may have im- 
plications in systems like the HTSC, with similar super- 
conducting parameters and possibly a similar mechanism 
for their superconductivity. It would be interesting to 
extend these measurements to other iron-pnictides with 
other type of substitution and doping levels, and also in 
larger samples to confirm whether the onset of the pre- 
cursor diamagnetism is close to 1.7T C , as proposed by the 
total-energy cutoff scenario. 
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